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Spatial variation of the fine-structure parameter and the 
cosmic microwave background 

Kris SigurdsonQ Andriy KurylovQ and Marc Kamionkowsk{| 

California Institute of Technology, Mail Code 130-33, Pasadena, CA 91125 

We study the effects on cosmic microwave background (CMB) temperature and polarization 
anisotropies of spatial fluctuations of the fine-structure parameter between causally disconnected 
regions of the Universe at the time of recombination. Analogous to weak gravitational lensing, 
in addition to modifying the mean power spectra and inducing a curl component (B mode) to the 
polarization, spatial fluctuations of the fine-structure parameter induce higher-order (non-Gaussian) 
temperature and polarization correlations in the CMB. We calculate these effects for the general case 
of arbitrary correlation between temperature fluctuations and fine-structure parameter fluctuations, 
and show the results for a model where these two types of fluctuations are uncorrelated. The 
formalism we present here may also be applied to other modifications of recombination physics 
that do not significantly alter the evolution of the dominant density perturbations. We discuss the 
constraints on the effective Lagrangian for variable fine-structure parameter necessary to realize this 
scenario. 

PACS numbers: 98.80.Es,98.70.Vc,06.20. Jr 



I. INTRODUCTION 

The possibility that the fine-structure parameter a 
may vary in time has long been entertained Q, 0, 0, 
IE S El and has received renewed interest with re- 
cent evidence from quasar spectra that may support a 
variation of less than one part in 10 4 over a time scale 
of ~ 10 Gyr HG3- Although the results may still be 
controversial, the observational work has inspired the- 
oretical work that investigates models with variable a 
mEEQHE[Hlll,as well as other work that 
investigates possible connections with dark energy and 
new long-range forces It has also stimulated 

a more careful re-investigation of the constraints placed 
on variable a from big-bang nucleosynthesis 22]. 

Of course, if a relativistic theory allows for temporal 
variation of a then it must also allow for spatial variations 
of a between regions not in causal contact. In this pa- 
per we study cosmological probes of spatial variations in 
a, focusing in particular on the cosmic microwave back- 
ground (CMB), which is rapidly becoming an increas- 
ingly precise probe of cosmological models [23, l24l l25j , 
as well as the physics that underlies them [2g. As we 
will sec, spatial variation in a induces a spatially-varying 
power spectrum. This induces "non-Gaussian" signa- 
tures in the CMB, in the form of locally anisotropic 
correlation functions, that cannot be described by the 
power spectrum alone (although strictly speaking the 
joint probability distribution of temperatures at n points 
remains a multivariate Gaussian). Quite interestingly, 
these effects are analogous to those of weak gravitational 
lensing ("cosmic shear") on the CMB, but differ in detail. 



So, for example, spatial variation of a can alter slightly 
the CMB power spectrum, induce a curl component (B 
mode) in the polarization, and induce higher-order tem- 
perature and polarization correlations in the CMB. Ef- 
fects like those we investigate here may also arise if there 
are other spatial variations in recombination physics that 
do not involve significant density/pressure perturbations. 
Our calculations are thus illustrative theoretically, apart 
from the specific application on which we focus. 

Below, we first review previous work [27H28l | that shows 
how the recombination history depends on a homoge- 
neous shift in a and then discuss how it affects the CMB 
power spectrum. We then show how spatial variations of 
a, in which the mean value of a is unaltered, affect the 
CMB temperature and polarization power spectra, and 
in so doing show that a curl component is induced in 
the CMB polarization. We then calculate the CMB bis- 
pectrum and trispectrum induced by spatial a variation. 
Throughout, we compare with the analogous calculations 
for weak lensing, and show how the effects of weak lens- 
ing and spatial a variations differ. We then discuss the 
properties of and constraints on a toy field-theory model 
for spatial a variation that produces the CMB effects 
we investigate here, without inducing significant density 
perturbations. 

Before proceeding further, we clarify that here we in- 
vestigate spatial variations in the fine-structure parame- 
ter a — e 2 /he that arise only from spatial variations in 
the electromagnetic gauge coupling e; we do not tinker 
with relativity nor quantum mechanics. 



II. RECOMBINATION AND a 
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Recombination depends on the value of a because the 
visibility function, the probability distribution of when a 
photon last scattered, is dependent on a. The visibility 
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function is defined as 



M -r dr 
9{t) = 6 dt 



where 
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is the differential optical depth of photons due to Thom- 
son scattering. Here, 
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is the Thomson cross section, n p is the total number den- 
sity of protons (both free and bound), and x e is the frac- 
tion of free electrons. The strongest effect of variations 
of a on this quantity occur due to the alteration of the 
ionization history x e (t). 

The recombination of hydrogen cannot proceed 
through direct recombination to the ground state because 
the emitted photon will immediately ionize a neighboring 
atom with high probability. Instead, the ionized frac- 
tion decreases mainly through the two-photon process 
2s — > Is, or via the cosmological redshifting of 2p — > Is 
Lyman-a photons out of the Lyman-a line. These pro- 
cesses are described by a single differential equation [2ij , 



dx e 
~dT 



j3e 



(1 — x e ) — 1Zn p x e 
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where /3 is the ionization coefficient, TZ is the recombi- 
nation coefficient, C is the Peebles efficiency factor (dis- 
cussed below), and 
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is the binding energy of the level with principle quantum 
number n. 

Through the Einstein relations we relate /3 to TZ 
through 



m e k B T 
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The recombination coefficient can be written 

oo n—1 
n=2 1=0 
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where 
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is the rate at which atoms recombine to the n, I energy 
level and w n is the efficiency for an n level to survive 
in a plasma |30j . The details of w n are not important 



for the present discussion, other than to note that at the 
densities of interest it is unity for n < n max ~ 500. Since 
the dominant contributions to 1Z come from n < 50, 1Z 
is insensitive to the weak a dependence of n max . Above, 
a n i is written in terms of the ionization cross section tr^{ , 
which can be expressed in the form [8ll | 



bf -ir (hv\ 
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Thus, we can write 



= HMf^]- e e/, 
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where 
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and it immediately follows that the a dependence of 1Z 
is of the form, 

K = a~ l TiF(£) = a 2 G(£,) . (12) 
As shown in Ref. |32j, for the temperatures of interest, 
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13V3. ( Si 
t>2 — —~ — In 



16tt 



kuT 



(13) 



(14) 



Given the scaling in Eq. (|12(l . we can then read off the a 
dependence. Explicitly, the recombination coefficient is 



K = 



3 y/2Trmlk B T \ 2k B T 



(15) 



The rate of recombination is inhibited by ionizing pho- 
tons which can disrupt atoms in the n = 2 state before 
they can decay to the ground state. The efficiency of 
recombination from the n — 2 state is described by the 
Peebles efficiency factor, 



C 



A H + A 2s ^i s 
A H + A 2s ^ u -+ 



(16) 



which is just the ratio of the recombination rates to the 
sum of the recombination and ionization rates from the 
n = 2 level. 
In Eq. (PJ), 



A H = 



8ttH 



(A 2p ^i s ) 3 n 1 



(17) 
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FIG. 1: The probability distribution for the last scattering of 
a photon, the visibility function, as a function of conformal 
time r\ in the ACDM model for tp = (a — ao)/ao = 0.03 
(dotted), tp — (solid), and tp — —0.03 (dashed). 
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FIG. 2: The angular power spectrum of CMB anisotropies 
for a spatially uniform tp — 0.03 (dotted), <p — (solid), and 
ip = -0.03 (dashed). 



is the rate at which a recombination is successful be- 
cause the emitted Lyman-a photon is redshifted out 
of the Lyman-a line before ionizing a hydrogen atom. 
Here, H = (l/a)(da/dt) is the Hubble expansion rate, 
m s ~ (1 — x e )n p is the number density of atoms in the 
Is state (almost all hydrogen atoms are in the Is state), 
and 

*2 P ^ls = 7, n 3 (18) 

is the Lyman-a rest wavelength. 

The two-photon process 2s — > Is proceeds through vir- 
tual atomic states at a rate A2 S _>i s = 8.22458 s _1 [33|. 
and scales as |34| 

A 2s ^ is oc a 8 . (19) 

Eqs. (|4"|>- l|19f) account for the a dependence of x e (t). 
Along with the a dependence of cr-r, this completely de- 
termines how g(t) varies with a in a given cosmology. 
While we expect a more complete calculation of recombi- 
nation, such as that in Ref . |35j , may yield further refine- 
ments to the a dependence, modifying hydrogen recom- 
bination as described above is adequate for our purposes 
because we are primarily concerned with the a depen- 
dence of the visibility function at very high redshift where 
the simple calculation accurately captures the physics. 
Also, it was determined in Ref. [28| that other effects, 



such as modifications to the details of helium recombi- 
nation, or the cooling of baryons, are small compared to 
the effect of variations of a on hydrogen recombination. 

In this paper we work within the flat geometry ACDM 
cosmology with baryon and matter densities Qb = 0.05 
and f2 m = 0.30, Hubble parameter h = 0.72, and spec- 
tral index n = 1. Fig. ^ shows the visibility function 
g[rf) = exp(— T)dr/dr] plotted versus the conformal time 
77 = J dt/a for three different values of tp = (a — ao)/ao 
where a = 0.00729735 ~ 1/137 is the value of the elec- 
tromagnetic fine structure parameter |36j| . For positive 
values of tp the visibility function is narrower and peaks 
earlier, while for negative values of ip the visibility func- 
tion is broader and peaks later. These effects impact 
directly the CMB angular power spectrum because the 
peak of the visibility function determines the physical 
distance to the last scattering surface, while the width 
of the visibility function determines the thickness of the 
last scattering surface. 

Fig. |21 shows the angular power spectrum of CMB 
anisotropies calculated assuming a spatially homoge- 
neous value of tp. The dependence of the spectrum on 
tp is easy to understand qualitatively. For positive tp the 
angular-diameter distance is larger, so the features are 
scaled systematically to higher values of I. The last- 
scattering surface is also narrower so that small-scale 
(high-/) features are less damped due to photon diffusion. 
For negative tp the opposite holds, the smaller angular- 
diameter distance scales features to lower values of I while 
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the broader last-scattering surface leads to more damp- 
ing of power on small scales. In the next two Sections we 
derive the changes in the angular power spectrum due to 
spatial fluctuations in <p between causally disconnected 
regions of the Universe. 



III. POWER SPECTRA 
A. CMB Power Spectra Fundamentals 

The CMB radiation is observed to be a nearly isotropic 
background of blackbody radiation at a temperature of 
Tomb = 2.728±0.004K |37|- Anisotropics in the temper- 
ature are observed with a fractional amplitude of ~ 1CP 5 
|38j |. and in the polarization with a fractional amplitude 
of ~ 10" 6 H^. For a review of the physics of CMB 
anisotropies see, e.g., Ref. |40j . 

The fundamental CMB anisotropy observables are the 
Stokes parameters (Q,Q,U,V), which can be expressed 
in the Pauli basis as (see, for example, Ref. [4l|) 

P(n) = 9(n)l + Q(h)cr 3 + U{n)a x + V(h)a 2 . (20) 

In this expression, 

T(n) - T CMB 



and 



0(n) 



TCMB 



(21) 



denotes the fractional temperature anisotropy in a direc- 
tion n, and the remaining Stokes parameters are normal- 
ized to this quantity. Here Q and U describe independent 
linear polarization states, while V describes circular po- 
larization. Because circular polarization cannot be gen- 
erated via Thomson scattering, V = for the CMB. 
It is convenient to introduce the quantities, 

± A(n)=Q(h)±iU(n), (22) 

which have the spin-2 transformation properties, 



z A(n) 



A(h) 



under a counterclockwise rotation of the coordinate axis 
by an angle <fi. 

We can expand 9(n) and ±A(n) in normal modes as 

El, 



7^( ± A M W)e ikX . (27) 



Here sY™ are the spin-s weighted spherical harmonics 
H3, with Y™ = oYf 1 , and A, (m) (fc) is the contribution 
to the angular mode Xi m from wave vectors of the pri- 
mordial density field of magnitude k. 

It is conventional to write the polarization in terms of 
the moments of the curl-free (scalar) configurations E\ m , 
and the moments of the divergence-free (pseudo-scalar) 
configurations Bi m , where 



-A, 



Eim. ± iBi. 



(28) 



We can then provide a complete description of an arbi- 
trary CMB anisotropy field using the moments 0/ m , Ei m , 
and Bi m . 

The basic observables of the random fields X(h) are 
the power spectra C? x , defined by, 



XX 



(29) 



where X, X £ {Q,E,B} and the angle brackets denote 
an average over all realizations. Here, 



C, 



xx 



tt(2Z + l) 2 



dk 



£ k 3 X^*(k)xl m \k) 



m=-2 



(30) 

A set of Gaussian random fields — and we expect 
{O, E, B} to be Gaussian — are completely characterized 
by their power spectra and cross-power spectra. Because 
the pseudo-scalar B has opposite parity to the scalars 
and E the only non- vanishing power spectra are Cf 9 , 
Cf E , Cf E , and Cf B . 

For small patches of sky it is an excellent approxi- 
mation to treat the sky as flat and expand the field 
X(h) in spin- weighted Fourier modes rather than spin- 
(23) weighted spherical harmonics. Thus we have (for exam- 

- . m 



pie Ref. |4J) 



6(n) 



d 2 \ 



9(l)e iln : 



(31) 
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and 



z A(h) 



where 
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Air 
21 + 1 



6 /m ^ m (n), (24) 
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(±A lm )( ±2 YT(n)) 



d 3 k 

(27j3" J 



e, (m) (fc)e ik ' x 



(25) 
(26) 



-A(h) 



r/ 2 l 



r±A(l)e 



±2i(0i-0) ihn 



(32) 



{2*r 

and we again define E and B through 

±A(\) =E{\) ± iB(l) , (33) 
In this notation the power spectra are defined by 

(A(l)A(l')) = (2^) 2 <5 2 (1 + l')Cf * • (34) 

Unless otherwise noted, we work within this flat-sky ap- 
proximation for the remainder of the paper. 
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B. Derivative Power Spectra 



How do the expressions for the power spectra change 
if we allow for spatial fluctuations of a? As a warmup we 
first consider a spatially uniform variation, a = ao(l+p), 
where ip <C 1. 

For a given primordial density field £(x), the tempera- 
ture and polarization patterns, 0(n) and ±A(h) 1 can be 
calculated by solving the combined Einstein equations 
and radiative-transfer equations, as well as the equations 
for the recombination history. As discussed above, this 
recombination history depends on a. Thus, the temper- 
ature and polarization fields are implicitly functions of 
<p = (a — a )/a . We can expand 0(n) = 6(n; <p) and 
±A(h) = ±A(h; ip) in Taylor series about ip = 0, 

6(n) = 9 (n) + d v e Q (h)cp + i^6 (n)^ 2 + ■ • • , (35) 



± A(n) = ± A (n) + d v (±A (n)) p + -d 2 (±A (n)) V 2 

+ ■■■ . (36) 
Inverting Eqs. I|31|) and (|32|l . we find that 

6(1) = J d 2 hO(h)e~ lhil , (37) 

and 



d z h{ ± A{h))e 



±2i( 



(38) 



Thus, the ip expansions can be written in 1 space as 

e(i) = e (i) + a v e (i) v + ^ 2 e (V + • • • , (39) 



-A(l) = ±A (l) + d v (±A (l))<p+-dl (±A (1)) p 2 + • ■ ■ , 



(40) 



(41) 



and in fact for any field X 6 {0, .B} we may write 
X(l) = X (l) + 5 v X (l)p + ^d 2 X (l)p 2 + • • 
To 0{ip 2 ) we then have, 
(X(l)X(l')) = (Xojl)Xo(l')) 

+ (x (i)^Xo(i'))p + (a v Xo(i)Xo(i'))p 
+ i(A (i)a 2 lo(i'))p 2 + i(a 2 x (i)x (i'))p 2 

+ (9 v Xo(l)^A (l'))p 2 , (42) 
and in terms of power spectra this becomes 



C xx = C x 



0X0 + (c? odx " + cf XoX "^ ip 

1 „y„;? 2 v„ 1 ^fi 2 : 



+ { ^Cf n9x " + ±Cf x ° x ° + cf XodXo l 



(43) 
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FIG. 3: The cf e ° se ° (solid) and cf° d2@0 (dashed) derivative 
power spectra. 

Since we are for the time being assuming a spatially 
uniform value of ip we may also write 

1 



xx 



d v C xx 



ip- 







2 ffr xx 
2 d * Cl 



(44) 



This allows us to make the identifications 



xx 



C, 



d v c xx 

d%C xx 



c, 



Xq Xq 



_ qXqOXq _|_ qOXqXq 



(jXqO Xq 



c l 



d 2 X X 



2(j^Xq8Xq 



(45) 



These identifications make it clear how to calculate the 
individual 'derivative' power spectra. Wc just differenti- 
ate Eq. (|30|l with respect to ip, evaluate the expression at 
(p. = 0, and pick off the terms with the requisite structure. 
For example, 

2 r dk 
~~k 



c: 



■dXgdXo 



7T(2Z + l) 2 



£ k 3 d v xi m) *(k)d v xr\k) 



m=-2 



(46) 



1,3=0 



where factors like dmX^' {k) can be calculated numer- 
ically or directly from first principles using the expres- 
sions for the x[ m \k) derived in, for example, Ref. |45j . 
Because they are used in subsequent calculations we 
have created a modified version of the code CMBFAST ^(J 
that can compute these derivative power spectra. (See 
Figs. 3-5.) 
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IV. SPATIAL VARIATIONS OF a 
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We now consider the effects of spatial variations of a 
(parametrized by (p) between different causally discon- 
nected regions of the Universe. 

First suppose that there were no density fluctuations, 
but spatial variations of a. In that case, photons from 
different points on the sky would be last scattered at dif- 
ferent cosmological times, but they would all still have 
the same frequency when observed by us. However, if 
there are density fluctuations, the manner in which they 
arc imprinted on the CMB depends on the value of a, 
as discussed above. Thus, if there are spatial variations 
in a, the power spectra (or two-point correlation func- 
tions) will vary from one place on the sky to another. 
This implies that the stochasticity of the spatial varia- 
tions in a induce non-Gaussianity in the CMB quantified 
by non-zero (connected) higher order correlation func- 
tions (trispectra and perhaps bispectra). It also implies 
a correction to the mean power spectrum (i.e., that mea- 
sured by mapping regions of the sky that contain many 
coherence regions of a), as well as the introduction of 
a non-zero curl in the polarization. All of these effects 
are analogous to similar effects induced by weak lensing 
of the CMB. The only difference is that in our case, the 
temperature and polarization patterns are modulated by 
a variable a, rather than lensing by an intervening den- 
sity field along the line of sight. 

In this Section, we first calculate the modified power 
spectra Cf e , Cf E , C EE , and Cf B . We then determine 
the form of the higher order correlations (bispectra and 
trispectra) in the next two Sections. 



A. Observable modes in the presence of <p 
fluctuations 

We assume that at a given position n at the surface 
of last scatter, the value of a is a(n) = a [l + <p(n)]. 
Here we treat ip(n) as a random field with angular power 
spectrum (p(l)p(V)) = (2tt) 2 5 2 (1 + l')Cf in the flat-sky 
approximation. 

We assume that the surface of last scatter is much thin- 
ner than the spatial correlation length of p, and that in a 
given direction a is constant throughout recombination. 
We also assume that the dynamics responsible for the 
variations of ip have a negligible effect on the perturba- 
tion evolution so that the sole effect of variations of p are 
a modification of the microphysics. We will discuss the 
2000 validity of these assumptions later. 

Again we expand our fields 0(n) = 0(n;<^) and 
±A(n) = ±A(h; ip) in a Taylor series about <p = as 



FIG. 5: The Cf e ° SE ° (solid), cf° a2j5 ° (dashed) and cf> s2e ° 
(dotted) derivative power spectra. 



6(n) =e (n) + d v e (nV(n) 

+ i^e (nV 2 (n), 



(47) 
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±A(n) = ±J 4 (n) + d v (±A Q (n)) ^(n 
1 

2°* 



+ ^d 2 ( ± A (h))^(h), (48) 



where y(n) is now a function of position. 

By taking the Fourier transform of 6(n) we find 

e(i) = e„(i)+[(0„e o )*d (1)+^ [(^e )*w] (i) (49) 



where 



(i) 



21/ 



(2^0 



r X(l')V(l-l') 



(50) 



is the convolution of two fields X and ip, and 



[X*ip*A] (1) = 



d 2 r 

(2^ 



[X«/>](1')A(1-1') (51) 



is the double convolution of three fields X, ip, and A. 

Similarly, taking linear combinations of the Fourier 
transforms of ±A(h) we find that 

E(l) = E Q (l) + [(d v E )* c tp] (1) - [{d v B )* B <p] (1) 



[(c££ )*cW] (1) 



(52) 



and 



B(l) = B„(l) + [(0 v B o )*c¥>] (1) + [(^£ )*s¥>] (1) 

+ 5 [(^KH (i) 



- [($,Eo)*.<p*<p] (1), 



(53) 



where 



[X* C <A] (1) = J cos (2^)X(1')V(1 - 1') , (54) 



and 



[** s ^] (1) 



21/ 



■sin(20 r )A-(I / )V(l-l'). (55) 



are the even- and odd-parity spin-2 weighted convolu- 
tions of X and ip respectively. 

Examining these expressions we find that a given 
mode X(l) receives corrections due to the combina- 
tion of modes {d™X {h), <p{h), <f{h), ■ ■ ■ , <p(K)} such 
that X)i=o = ^ Furthermore, the E and 
_B modes mix so that, for example, the mode 
-B(l) can be induced by the combinations of modes 
{d^Eo%), <p(h), <p(h), ¥>(!*)} such that J2t h = I- 
These effects modify the angular power spectra of CMB 
anisotropics, and introduce higher-order connected (non- 
Gaussian) correlation functions. 



B. The 66 Power Spectrum 

Using Eq. (|49|l we find that the expansion for the two 
point correlation function (in Fourier space) is 

(6(1)6(1')) = (6o(l)6 (l')> 

+ (6 (1) [d v QoM (l')> + (K&oM (l)6o(l')> 

+ i<e (i) [a*e *^] (!')) 



-{[& v Qv*P*p] 0)00 (l')> 
([d v Q *<p] (1) [d v Q *p] (!')> 



(56) 



to 0(tp 2 ). In the above expression and what follows we 
adopt the convention that the differential operators d v 
act only the the field immediately following them. 

We assume that 60 and tp are zero-mean Gaussian ran- 
dom fields without higher-order connected correlators. 
By writing out the convolutions and Wick expanding the 
correlators it is easy to verify that correlators involving 
an odd number of fields vanish, and so there are no cor- 
rections to first order in tp. It is also straightforward to 
verify that 



(e (i)[^e ^(iO) = ([€6o*w](i)e (i')) 



= (2tt) 2 J 2 (1 + 1') 



where 



r (V-A) _ 



d 2 l 



(27T) 



G 



2^1 



(57) 



(58) 



is the covariance between two fields ip and A. Similarly 
we can show that 

([d v e *ip] (i) [d v B *<p] (!')) 

= (2tt) 2 5 2 (1 + l'){ [C^*C d6ode °] l 

+ [C deov *C 90ov ] l }, (59) 

where we have dropped terms that contribute only when 
1 = 0. 

Collecting all terms we find that to leading order the 
average power spectrum including fluctuations in ip is 

Cf 6 = C &0@0 + cj(w)c 0od2e " + 2a (vd2Bo) Cf°' f 

+ [C^-kC d@ad@0 ] l + [C deoip -kC 9e ° v ] r 

(60) 

Note that the corrections to involve couplings be- 

tween the derivative power spectra, which are calculated 
as described in Section UlI Bl and Cf v and the various 
cross power spectra, which are specified by the model 
that generates spatial variations in tp. 

In general 60 and tp may be correlated if, for instance, 
they are generated by a common mechanism or if they 



8 



5x10-1° I- 




In this case the average 00 power spectrum is 

■(<p<p) 



2000 



(65) 



where /„ is the n th -order modified Bessel function of 
the first kind. We show this average power spectrum 
in Fig. |SJ The main effect of ip fluctuations is to reduce 
the amplitude of oscillatory features in the damping tail. 
This effect can be understood by noting that patches of 
the sky with different values of ip will have different power 
spectra, and that, as can be seen in Fig. [3 the location 
of of the peaks in the damping tail of these power spectra 
shift as ip changes. These patch power spectra add inco- 
herently, so that the amplitude of the oscillatory compo- 
nent of the average power spectrum is reduced. This is 
the same type of effect as in weak gravitational lensing 
|47|. As we will see, this is the predominant effect in the 
other power spectra as well. 



FIG. 6: The solid curve shows the OO power spectrum for the 
Gaussian correlation function for tp with angular correlation 
scale 6 C = 1° and variance = 9 x 10~ 4 (fluctuations in ip 

at the 3% level). The dashed curve shows the power spectrum 
without fluctuations in ip. 



are strongly coupled through evolution equations. We 
discuss in Section IVTT1 why we do not expect the latter 
source of correlations to be important as long as the en- 
ergy in the ip field and its fluctuations are small compared 
to the dominant radiation and matter perturbations. In 
the case where ©o and ip have no cross-correlation, the 
expression simplifies to 



ee 



c; 



e e 
i 



9e ae i 



(61) 



The result for Cf will of course depend on Cf v . To 
illustrate we consider a simple model in which ip is highly 
correlated on angular scales smaller than the correlation 
angle C , and uncorrelated on larger scales. (Below we 
discuss a physical model that may produce such a corre- 
lation function.) We thus have, 



{<p(0)<p(P)) =o -(w) e -(<V0c) 2 
In the flat-sky approximation we have 



Qipifl _ 



which implies that 



d 2 e (<p(pM0)) 



-no 



(62) 



(63) 



(64) 



C. The QE Power Spectrum 

Using Eqs. (|4*9"|> and (|52*)l . and immediately dropping 
the vanishing correlators involving an odd number of 
fields, we find that the expansion for the power spectrum 



(e(W)) = (e (iW)> 
i 



(0 O (1) [fl£Eb*c*H (!')) 



-<e (i) [^H (10) 
<[^e ^] (i) [a v Eo*M (10) 

([d v QoM (1) P„flo*.d (10) 

i([^e *w] (i)^o(iO) 



(66) 



to 0{<p 2 ). 

By writing out the convolutions, Wick expanding the 
correlators, and noting that terms involving Bq vanish 
due to parity, we find that the non-vanishing terms are 



(6 (1) [fcEoMpup] (10) 



(2^)^ 2 (1 + 1' 



(67) 



[c£e *¥H (i)£b(iO) 

= (2tt) 2 J 2 (1 + 1') 



(68) 
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(72) 



We show this average power spectrum in Fig [7] The 
major effect of tp fluctuations on the QE power spectrum 
is to reduce the peak amplitudes on small scales. 



D. The EE Power Spectrum 

Using Eq. (|52[1 , and dropping the correlators involving 
an odd number of fields and those that vanish due to 
parity, we find that the expansion for the power spectrum 



is 



FIG. 7: The QE power spectra for the Gaussian correlation 
function for tp with a {tptp) = 9 x 10~ 4 and 6 C = 1° (solid) and 
9 C — 2° (dotted). The dashed line shows the power spectrum 
without fluctuations in p. 



([d v Qa*<p] (1) [9 V E *M (1')) 

= {2n) 2 5 2 (l + l'){ [C ae ° dEo * c C w ] l 
+ [C a£ov * c C ae ° v ] ; }. 



(69) 



Collecting these terms we find that, to leading order, 
the average cross power-spectrum including fluctuations 
in tp is 



qQE _ (jQoE _|_ ^(tptp) 



1 r e a 2 E 1 r E d 2 e 



+ (J (vd 2 E ) ( j<3> <p + (J {<pd' 2 Q ) ( jE Q9 

+ [c de ° dEo * c c vv ] l + [c 9Eov * c c de ° v ] l . 

(70) 



If tp has no correlation with the density field, the ex- 
pression simplifies to 



(E(l)E(V)) = (£o(l)£o(l')> 

+ ~<£ (1) [d 2 v E * c ^] (1')) 

+ \{[d 2 v E^ c ^] (1)£ (1')) 
+ ([d v Eo* e <p] (1) [d v Eo* e <p] (1')) 
+ {[d v B * B <p] (1) [d v B * a <p] (l')> 



(73) 



to 0{ip 2 ). 

After evaluating the convolutions and Wick expanding 
these correlators we find 



(Eo{l) [d 2 p E * c tp*tp] (l')> = ([^3>*«W] {\')E Q {\)) 
= (27r) 2 (5 2 (l + l' 



a ( w ) (jEod 2 E j l _2 a (<Pd 2 E )QEoip 



(74) 



[d v Eo* e <p] (1) [d v Eo* B <p] (1')) 

= (2n) 2 5 2 (l + Y) f ^[ C os 2 (2^)Cf^C^ 



+ cos (2<f> h ] 



2lSsm z (0|J 
ll-lil 2 



- 1 



s-ydE tp ^tdEo tp \ 



(75) 



(J®E _ qQ Eq _j_ a (tptp) 



1 



de dE ^ 



(71) 



Inserting the power spectrum from Eq. (|64|l , we obtain 



( [d v B * s ip] (1) [9 VJ B *s<d (l')> 
= (2.) 2 <5 2 (1 + 1') J ^ s in 2 (2^Cf^C^ 



(76) 
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FIG. 8: The EE power spectra for the Gaussian correlation 
function for tp with a {tptp) = 9 x 1(T 4 and 6 C = 1° (solid) and 
9 C — 2° (dotted). The dashed line shows the power spectrum 
without fluctuations in tp. 



Collecting these terms we find that 



C-, 



EE 



qE E _|_ a (tpip)(jE d 2 E _|_ 2 a (<pd 2 E ) qE oV 



'I 
d 2 V 



(2 
2 ( 

cos (2(j) V ) 



-L{[cos 2 (2^)C? E ° dE ° 

7T) <• 



+ sin 2 {2^)Cj > , BodBo ]C^ 1 , l 



2r 2 sin 2 (<j> v ) 
ll-l'l 2 



- 1 



c dE oVc dEo^ ^ (77) 



is the expression for the average 5-mode power spectrum 
including fluctuations in tp. 

If Eq and iy9 are uncorrelatcd, the expression simplifies 

to 

(jEE (jEoEo _|_ g-iw) (J^o d 2 Eq 

d2V r cos 2 (2^)^ B ° aB ° 



(2tt 



+ sin 2 {2(f> l >)Cf, B ° dBo }C^ i , l 



(78) 



For the tp power spectrum in Eq (|64|l . dropping the 
negligible primordial 5-mode term, we obtain the ex- 



C EE = C E E + a{w ) J C Eod*E + ^ j tf^-iV+l")'* 



2 



4 ( f U 



9 2 



(79) 



We plot this average power spectrum in Fig. |H1 Again, 
the effect of tp fluctuations is to reduce the amplitude of 
the oscillatory component at small angular scales. 



E. The BB Power Spectrum 

Using Eq. Ij53(l . and once again dropping the correla- 
tors involving an odd number of fields and those that 
vanish due to parity, we find that the expansion for the 
two point correlation function is 

(5(1)5(1')} = (50(1)50(1')) 

+ \{B Q {\) [S 2 i?o*c¥H (!')) 

+ ~{[dlB * cV *tp] (1)50(1')) 

+ ([d v B * c tp] (1) [d v B Q * c tp] (1')) 

+ {[d v E -k s <p] (1) [d v E * s tp] (1')) (80) 

to 0{ip 2 ). 

Evaluating the convolutions and Wick expanding these 
correlators leads to the expressions 

(5 (1) [9 2 5 * c ^] (1')) = ([d 2 5 * c ^] (l')So(l)) 
= (2tt) 2 (5 2 (1 + 1') L<-v<p) C Zod 2 B ] ^ ( 81 ) 



; [8 V B *M (1) [d v B * c tp] {!')) 

= (2.) 2 5 2 (1 + 1') / |3L cos 2 (2<P h )C? B ° SB °CW li{ 



(82) 



[d v Ea* B (p] (1) [d v Ea* m tp] (I')) 

= <2k)'P(L + Y) / |^{siu 2 {2<t> h )Cf B ° dE °C^ 
+ sin(2^ 1 )^ ?Sin(2 ^ ) - MlSin( ^ ) " 



dE Q ip s-idE: 



|l-li| 

'ov\ 



(83) 



Any </?-5o cross correlations must vanish due to parity. 
Collecting all terms we find that the expression for the 
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FIG. 9: The BB power spectra for the Gaussian correlation 
function for if with a {{fap) = 9 x 1(T 4 and C = 1° (solid), C = 
2° (dotted), 6 C = 5° (dashed), and 9 C = 10° (dot-dashed). 



average _B-mode power spectrum is 

jdB 



d 2 V 



{2-Kf 



{[cos 2 {2^)Cf, Boi 



+ sin 2 (2<f> ll )Cf l EodEo 



sin (2^/) 



r«p<p 



I' 2 sin(20j/) - 2/i' sin I 



1-1' 



(84) 



and if £o and (/? are uncorrelated the expression simplifies 
to 

qBB _ qB B _j_ a ( W )(jB a d 2 Bo 

d 2 V 



d B d B 



(2nf 



[cos' {2<j> v )C?, 
+ sin 2 (20,/)C^ BoaBo ]^„ 



i-i'| ■ 

(85) 



If the intrinsic B-modes due to gravitational waves are 
negligible then for our tp power spectrum, Eq. (|64(l . the 
average B-mode power spectrum is 



C BB = \ I d l'l'e'^ l2+l ' 2 ^ 2 



2 

c ;;/ 



h [ fir 



QdE dEo | 



(86) 



This expression shows that E modes modulated by ip 
fluctuations can induce B modes. We note here that this 
effect is more general than the specific application of a 
variation we focus on here, and is a generic feature of 
any fluctuations that modulate the power spectrum but 
have negligible effects on the evolution of the dominant 
density perturbations. In Fig. 03 we plot the induced B- 
mode power spectra due to ip fluctuations. The induced 
-B-mode power spectrum inherits the oscillatory features 
of the unperturbed E'o-mode power spectrum as long as 
correlation angle is larger than the horizon size at re- 
combination, while the amplitude of the power spectrum 
decreases as the correlation angle increases. 



V. BISPECTRA 

We have shown in the previous Section that modula- 
tion of the temperature and polarization power spectra 
by spatial variations of a alters the mean power spectra 
of CMB anisotropies. In this Section and the one that 
follows we show how this modulation introduces higher 
order correlations into the CMB temperature field. 

The temperature bispectrum is defined in terms of the 
connected piece (the terms remaining after the Gaussian 
piece is subtracted out) of the three point correlation 
function in Fourier space as 



(e(i 1 )e(i 2 )e(i 3 )) ( 



= (27r)^(li + h + l 3 )S www (li, la, Is) • (87) 

This expression must be invariant under the exchange 
of any two fields, or equivalently, 1 vectors. We insert 
Eq. (|49|l into the three-point correlation function, and 
after some straightforward algebra we find that 



-B° ee (li, I2, 13) — ^ B(li,lj,l6_j-j0 

i j=l 



where the sums runs over the distinct permutations of 
{1, 2, 3} and, to leading order, 



Bft.lj.lfc) = cf a9e °Cf 0V + cf olp Cf ode ° 



(89) 



Thus the temperature bispectrum vanishes unless tp 
and Oo are correlated. In the models we consider in this 
work, tp and Oo are not expected to be highly correlated. 
However, if a model of spatial variations of a predicts 
that variations of tp and Oo are strongly correlated the 
bispectrum may be a strong signature of such a model as, 
unlike the power spectrum or the trispectrum (discussed 
below), it is first order in tp. We note that expressions 
analogous to Eqs. I|88|) and l|89|) will hold for the polar- 
ization and cross-bispectra as well. 
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VI. TRISPECTRA 

In analogy with the bispectrum, the temperature 
trispectrum is defined in terms of the connected piece 
of the four point correlation function in Fourier space as 

(e(ii)e(i 2 )e(i 3 )e(i 4 )) c 

= (2TT) 2 S 2 (h + h + h + l 4 )T eeee (l 1 , I a , la, 1 4 ) . 

(90) 

This expression must also be invariant under the ex- 
change of any two fields, or equivalently, 1 vectors. We 
insert Eq. I|49|l into the four-point correlation function, 
and after some straightforward algebra we find that 

4 4 

T eeee (li) l2) h u) = £ £ ^(1,, 1„ l fc , L) 

i<3 k<l 
4 4 

+ J2 J2 TB(ii,i,-,ife,io, 



i,j,k=l 1=1 
i<j<k l^i,3,k 



(91) 



where the sums runs over the distinct permutations of 
{1,2,3,4}, 



*-A{k, ij, lfc, H) — L>i } ^ll^+lfcl + 

_l_ / ^>0o90o^f0oV I (jQo<P (j&od&a\ ( qSOop ^ £id& n <f 

\ U l j U h J \ I lj +lfc I 

+ H C h \ C \h+h\ ) ' V* 2 ) 

and 

T B (l<,l;,Ik,Ij) = 2C° od20o C° olp C° oV + 

+ 2C°"' p (C° oip C° od200 + C® od20o C° oip ) . (93) 

If Oq and ip have no cross correlation, the second term 
vanishes and the first term simplifies to 

rp /i l l 1 \ _ ,o0o 90o/^0o 9©o I r<f"P \f~tff \ 

(94) 

Expressions similar to Eqs. H91f) - i|93|) wm hold for the 
polarization and cross-trispectra as well. 



A. The Kurtosis 

The trispectrum is a nontrivial function of 6 variables 
and extracting the full trispectrum from CMB data will 
be a challenging experimental endeavor. It is therefore 
worthwhile to examine simpler statistical signatures of 
non-Gaussianity in the CMB. If the trispectrum is non- 
vanishing the probability distribution function of 0(n) 
on the sky will no longer be precisely Gaussian. As in 



the case of weak lensing, the deviation from Gaussian- 
ity can be parametrized by the kurtosis [48|. which is a 
measure of how flattened out or peaked the distribution 
is relative to a Gaussian distribution. 

The kurtosis of a non-Gaussian random field may be 
written in terms of the trispectrum as 



K{6) 



1 f d 2 h d 2 l 2 d 2 h d 2 l 4 



a 4 (9) J (2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 
x {(2tt) 2 5 2 (1i + h + h + U)T eee@ (h, h, I 3 , h) 

x w{he)w{i 2 e)w{h6)w{U6)} , (95) 

where W(19) is a smoothing function with smoothing 
scale 9, and 



* 2 (0) = J ^C? e W 2 {W) 



(96) 



is the smoothed variance. 

If we adopt a Gaussian smoothing function 

W{W) = e-2<^ 2 



(97) 



with <jb — 0/(v81n2) and insert the expression for the 
trispectrum from Eq. (|94|l we find after some algebra that 

K{9) = ^4(0) / hdhl 2 dhhdl 3 {c® de C^C? 3 d@ 

x I {o%hh)Io(o%hh)e-'*tt + % + ®] } • (98) 

It was shown in Ref. that the sample variance of 
the kurtosis of a Gaussian random field smoothed over 
an angular scale 9 for a full-sky experiment is 



>K(8) 



(99) 



In Fig.EHwe plot the ratio of the kurtosis to the standard 
deviation of the kurtosis as a function of smoothing scale. 
Despite the improved variance at high resolution the kur- 
tosis due to weak lensing cannot be detected because in 
the limit of infinite resolution the kurtosis of weak lensing 
vanishes more quickly than the variance. This can be un- 
derstood by noting that weak lensing is power conserving 
because it maps the temperature at one point on the sky 
to another, and so with infinite resolution the probability 
distribution function is Gaussian. In contrast, the kur- 
tosis due to spatial a fluctuations approaches a constant 
value in the limit of infinite resolution, and so for low 
9 the signal-to-noise increases like 0~ 1 . This occurs be- 
cause a fluctuations do not lead to a pure remapping of 
of the temperature pattern (they do not conserve power), 
and thus in the limit of infinite resolution modulate the 
variance of the (Gaussian) temperature probability dis- 
tribution function from one patch of the sky to another. 
The resulting mean probability distribution is no longer 
Gaussian. Note that the dip in the kurtosis at 9 ss 15 ar- 
cminutes occurs because when smoothed over that scale 
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FIG. 10: The ratio of the kurtosis to the variance in the esti- 
mator of the kurtosis due to Gaussian fluctuations for spatial 
a fluctuations with a {vv) = 9 x 10~ 4 and 9 C = 1° (solid) 
and the same ratio for weak lensing (dashed). The kurtosis 
is undetectable for weak lensing at any angular resolution, 
while for spatial a fluctuations of this amplitude the kurto- 
sis may in principle be detected by a high resolution no-noise 
experiment. 



modulation of the CMB due to spatial a fluctuations be- 
comes approximately power conserving. An ideal noise- 
free experiment at an angular resolution of 1 arcminute 
could detect the kurtosis due to a fluctuations at the level 
of a^) = 9 x 10 4 , while higher resolution experiments 
could observe lower amplitude a fluctuations. 



B. A Discriminating Filter 

As discussed above, weak lensing of the CMB by mat- 
ter along the line of sight must induce a contribution to 
the trispectrum 



L eeee (i!, i 2 , i 3 , i 4 ) = E E L &> l s> 1 *. »o . 

(100) 



i<3 k<l 



where 

L(li,lj, lfe, 1;) 



£iOo©o£i6o©o ^ 
Ik h 



c | tf +lfcl [(i l + i fe )-i fe ][(i J + i fc )-i ; ] 



cj; +ul [(i J -+i fc )-i fc ][(i J -+i fc )-i 1 



(101) 



We emphasize that in this expression Ct is the power 
spectrum of the projected lensing potential and not the 
power spectrum of a fluctuations Cf v . 

Since this contribution to the trispectrum must be in 
the CMB, we now derive the a filter to distinguish be- 
tween the weak-lensing and a-fluctuation trispectra. 

For compactness we introduce the notation 



Xijk = X 



eeee 



(1, 



7 lj 7 Ifc 7 



L - I," 



and 



V= ( — 



d 2 h 



) 2 (2tt) 2 (2tt) 2 



(102) 



(103) 



The signal for detection of the a-fluctuation trispec- 
trum can be in general written as a windowed inte- 
gral over trispectrum configurations (quadrilaterals in 1- 
space) , 



s = y w t jkTij k 



ijk 



(104) 



where Wijk is the to-be-determined window function. 

To determine the least-square quartic discriminator we 
treat the trispectrum due to weak gravitational lensing 
as a source of noise that should be minimized and write 



N L = ^2w ilk L i:jk . 

ijk 



(105) 



In addition we at minimum also have the Gaussian 
noise due to cosmic variance 



ijk 



where 



Gijk 



4! 



4tt/ ; 



sky 



U U u h u h u \h 



and 



ee 



„ee , ^fskys 2 2j2 

l ex-p ± CMB 



(106) 



(107) 



(108) 



is the sum of the actual CMB power spectrum and the 
noise power spectrum introduced by an experiment that 
observes a fraction f sky of the sky with beam width 
<7b for a time t exp with detectors of noise-equivalent- 
temperature s. 

The total noise is then just 



N 



N% + N? 



(109) 



We want to solve for the W^k that maximizes the signal 
to noise ratio S/N. To do this we set 



5W ilk 



N 2 



0. 



(110) 
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After some straightforward algebra, and making use of 
the fact that S 2 /N 2 is invariant under renormalization of 
Wijk, we find that the signal-to-noise ratio is extremized 
only if Wijk takes the form 



Tijk — XL 



ijk 



G 



ijk 



Thus we have 
S 2 



T 2 - 2XTX + X 2 X 2 



N 2 T + X 2 - {2XX- X 2 C){l + C) 



where 



T = 



\ " 1 ijk 



X 



... Gijk 
ijk J 



ijk 



G 



T 2 

L ijk 

Gijk 



and the optimum weighting is 

X 



A = 



l + C 



(111) 



(112) 



(113) 



(114) 



(115) 



(116) 



In the absence of lensing this expression reduces to the 
conventional signal-to-noise measure 



c2 t 2 

_ rj. _ ^ ijk 



ijk 



G 



ijk 



(117) 



In Fig. II II we show the maximum signal-to- noise ratio 
for detection of the a-fluctuation trispectrum that would 
be detected were weak lensing not present. For a given 
weak lensing power spectrum Ct the actual signal-to- 
noise will be less than the bound shown. 



VII. THEORETICAL MODELS OF VARIABLE Oi 

Recent theoretical work [HI E3, E3, EE EE EE 

I20L I2H has considered the ingredients required of field- 
theory models for variable fine-structure parameter in or- 
der to explain a small time variation of a. Here we briefly 
amend those discussions to consider the field-theory re- 
quirements for spatial variations of the sort we consider 
in this paper. 

The simplest way to introduce spatial variation of a 
is to couple the photon to a scalar field <f>{x) through a 
term (see, e.g., Ref. [2(J and references therein) 
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FIG. 11: The signal-to-noise ratio of the a-fluctuation trispec- 
trum that would be observed by an experiment with s = 
12.42 pKy'sec observing for 1 year (the effective parameters 
of the Planck satellite) with resolution 6 if weak lensing were 
not present for 0-( vv ) = 9 x 10 -4 . Since in practice we know 
weak lensing must be present in the CMB this curve serves 
as an upper bound on the observable signal-to-noise ratio. 



in the Lagrangian that appears in addition to the usual 
electromagnetic Lagrangian, F^ v 'F^ v j '4. Here, is a 
constant with dimensions of mass, A is an overall cou- 
pling constant, and g(x) — X)nLi c nX n is a dimensionless 
function of the ratio 4>/ normalized so that <?(1) = 1. 
This interaction then leads to a fractional change in a, 



S a 
a 



-Xg 



(119) 



The best present constraint to g comes from the evolution 
of globular-cluster stars; this requires [X(dg / dcf))]^ 1 > 
1.6 x 10 10 GeV ^E|. If we consider fractional fluctuations 
in a of < 1CT 2 , then A - 1CT 2 , 6<j) < fy, and - A 2 . 

Cosmological mechanisms that might induce spatial 
variations in (f> are analogous to those that have been con- 
sidered, e.g., for spatial variations in the axion field [50| . 
These variations in <j) can arise either during or after infla- 
tion. If cj) is a spectator field during inflation, then fluctu- 
ations in cj) can be induced quantum mechanically during 
inflation resulting in a nearly scale-invariant spectrum for 
4> fluctuations. Alternatively, if cj) is a pseudo-Nambu- 
Goldstone field in a model with an approximate global 
symmetry, it could fall during the symmetry-breaking 
phase transition to different random points on the vac- 
uum manifold (which is periodic in <f> with period 2/^) 
in different causally-disconnected regions. In this case, 
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the scalar-photon interaction will most generally be an 
arbitrary Fourier series in <f>. A simple interaction of this 
form, containing only the first harmonic and (as moti- 
vated below) constraining g(<f>/ f<f>) to be an even func- 
tion, is 



C 



1 — cos — - 



F^F, 



[lis ■ 



(120) 



For such an interaction (f> will be uncorrelated on scales 
larger than the horizon, and da /a is fixed by the magni- 
tude of the explicit symmetry breaking term A, rather 
than by fj,. In either case (inflation or spontaneous 
symmetry breaking), the gradient term in the scalar- 
field Lagrangian will tend to align the scalar field within 
causally-connected regions of the Universe. Thus, the 
value of a at the surface of last scatter should be con- 
stant within square-degree patches, but will vary from 
one square-degree patch to another, as we have assumed 
throughout this paper. 

Now consider constraints to the potential-energy den- 
sity V(4>) for the scalar field. If it is quadratic, V(<fi) = 
m^0 2 /2, then (neglecting for the moment the interaction 
with photons) the equation of motion for <j) is 



<p + 3H(t) + m\4> = 



and has the solution 



<Kt) = 



sin(m0t) 



(121) 



(122) 



in the matter dominated era. At early times when 
H 3> m<j, the value of (f> is frozen at <f> w cf>Q, while for 
H < the field in each horizon will oscillate about 
its minimum with an amplitude decaying like 1/t. We 
thus require < H lcc ~ 10 -28 eV so that the a 
fluctuations are frozen in at recombination. On the 
other hand, observations of the Lyman-a forest at red- 
shifts z ~ 4 show no evidence of spatial fluctuations 
in a larger than one part in 10 4 . Thus, if we consider 
(5a/ a) > 10~ 4 at the surface of last scatter, then we 
must require that Xg\ z ~4 < 10~ 4 . If the leading order 
term in g is quadratic in <f> (as we argue below) this 
leads to the constraint > y/X/0.01 x 10~ 31 eV. Fi- 
nally, laboratory experiments constrain the value a/ a 
today to be < 5 x 10~ 15 yr _1 To satisfy this re- 

quirement we require that g cx <fi 2 to leading order at 
small <fi so that a/a oc t~ 2 . This leads to the constraint 
> 2 x (A/0.01) x 10" 30 eV if the oscillatory factor 
is unity. We note that if g oc 4> n with n > 2 for small 
(f> this constraint can be relaxed, and that late-time A 
domination does not significantly alter the limit. These 
constraints will insure that the model conforms to upper 
limits to a variations in the low-redshift Universe, but 
we caution that they are constraints based on the root- 
mean-squared value of the scalar field. If we happen to 
live in a region where the amplitude of </> is randomly 
small (large) the constraints based on local observations 



will be weaker (stronger) than discussed above. Also, the 
nonlinear evolution of the mass distribution may modify 
constraints based on recent terrestrial phenomena ■ 

We must also be sure that fluctuations in the scalar 
field do not lead to density perturbations that exceed 
those of amplitude 10~ 5 . This constraint requires that 
the gradient-energy density, k 2 (A(f>) 2 < 10~ 5 /O m (z rcc ), 
where k ~ H vcc is the largest wavenumber for which per- 
turbations are significant at the redshift z ICC ~ 1100 of 
decoupling, and iJ roc = fill 2 H Q (l + z ICC ) 3 / 2 is the Hubble 
parameter at decoupling. However, from the Friedmann 
equation, p m (z rcc ) ~ H 2 cc mp l7 where mpi ~ 10 19 GeV 
is the Planck mass, so we find a constraint A(f> < f$ < 
3 x 10 16 GeV. 

Now lets examine the effect of the scalar-photon La- 
grangian on the scalar field dynamics. If we include the 
interaction in the Lagrangian it is easy to verify the equa- 
tion of motion for the scalar field is modified to include 
a forcing term 



-^i. (p F» v ) - 



where 



(123) 



(124) 



and {E 2 ) and (-B 2 ) are the spatially averaged squared 
electric and magnetic fields respectively. A bath of ther- 
mal photons has {E 2X j = {B 2 ) and so will not contribute 
to the forcing term. However, nonrelativistic matter is a 
source for electromagnetic fields with (-E 2 ) ^= (B 2 ), and 
so the forcing term should be proportional to the density 
of nonrelativistic matter p m . If, as we discussed above, 
g(4>/ f(j>) is dominated by its quadratic term then in the 
matter dominated era the equation of motion reads 



4>+j<j)+(m 2 



0. 



(125) 



where r\ — (2A£p mo to)//J, and £ is the fraction of mat- 
ter density p m due to electromagnetic energy. 1 Thus, if 
X) is small the effect of the forcing term is to introduce 
a small time-dependent mass term into the equation of 
motion. The solution will qualitatively behave like that 
of Eq. I|122(l , except that the time that field begins to os- 
cillate may shift by a small 77-dependent factor, and the 
power law of decay will shift from —1 to — 1 + r\. Gener- 
ically r\ need not be small, but for > 5 x 10 15 GeV 
and A « 0.01 we find |7y| < 0.1 2 This brings us uncom- 
fortably close to the limits on from gradient energy 



1 Naively, one might guess that ~nG intergalactic magnetic fields 
would suggest g ~ 10 -11 . However, the quantity of interest here 
is (B 2 ), rather than (B) 2 , as well as (E 2 ). Thus, strong mi- 
croscopic magnetic fields in the vicinity of electrons and protons 
may contribute £5 ~ 10~ 5 |2(j . 

2 This assumes a baryon fraction Qf,/Q m = 1/6, that magnetic 
energy contributes 2 X 10 — 5 |2flj . and that the dark matter 
has negligible electromagnetic energy so that § = (£lb/£lm)£b- 



16 



density. In a forthcoming publication |53j we argue that, 
in general, rj should be replaced by the quantity fj = rj+uj, 
and that fj may be small even if 77 is not. The term ui 
arises from the terms of the effective non-renormalizable 
Lagrangian of the form X n m^,(j) 2n, ipip or <7 ra m 2 <A 2ra x 2 that 
couple the scalar field <p to the energy density in other 
matter fields. We note in passing that for those theo- 
ries that do not lead to 'fifth forces' which violate the 
weak-equivalence-principle (WEP) fj — precisely. 

We thus have a constraint [X(dg / 'dcj))}^ 1 > 10 10 GcV to 
the scalar-photon coupling, and a constraint (2A/0.01) x 
1CU 30 eV < rritf, < 10 -28 eV to the scalar-particle mass. 
We note that these relations may seem hard to recon- 
cile, as contributions to the scalar-particle propagator 
from divergent loop diagrams containing photons should 
generically be large. Likewise, if is a Goldstone mode 
for some spontaneously broken global symmetry, it may 
be unusual from the point of view of Planck-scale physics 
54]. It was also pointed out in Ref. [5^ that variations 
of a light scalar should give rise to large variations in vac- 
uum energy density. On the other hand, similar problems 
arise in models for inflation and for scalar-field models 
for dark energy and dark matter. Here, we take the view 
that our ^-photon coupling is simply a low-energy effec- 
tive Lagrangian, and anticipate that the light mass can 
be protected, e.g., by a mechanism such as that proposed 
in Ref. (5(| where the scale of an inflaton potential is 
fixed by the size of an extra dimension. If so then one 
might assume that whatever mechanism alleviates the 
cosmological-constant problem may also alleviate prob- 
lems associated with vacuum-energy gradients. 



time of recombination. Although we have focused on the 
particular case of fluctuations in the fine-structure pa- 
rameter, the formalism we have presented may be applied 
mutatis mutandis to other modifications of recombination 
physics that do not alter the evolution of the dominant 
density perturbations. As discussed above, such fluctua- 
tions will alter the predicted CMB power spectra, intro- 
duce a B-mode polarization signal, and introduce tem- 
perature and polarization trispectra and perhaps bispec- 
tra. We stress here that these results are not dependent 
on the particular model for a variation discussed in Sec- 
tion Ell 

From the point of view of effective field theory, varia- 
tions of the fine-structure parameter can be phrased in 
terms of a scalar-photon interaction Lagrangian >C^ 77 , 
the parameters of which can be chosen to be consistent 
with current experimental limits. Ultimately, if such 
light, cosmologically interesting, scalars do exist in na- 
ture their mass must be protected by some yet unknown 
mechanism. 
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VIII. CONCLUSION 



In this paper we have studied the effects on the CMB 
of spatial fluctuations in the fine-structure parameter be- 
tween causally disconnected regions of the universe at the 
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